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Abstract 

Hans Bethe contributed in many ways to our understanding of the supernovae that 
happen in massive stars, but, to this day, a first principles model of how the explo- 
sion is energized is lacking. Nevertheless, a quantitative theory of nucleosynthesis is 
possible. We present a survey of the nucleosynthesis that occurs in 32 stars of solar 
metallicity in the mass range 12 to 120 Mq. The most recent set of solar abundances, 
opacities, mass loss rates, and current estimates of nuclear reaction rates are em- 
ployed. Restrictions on the mass cut and explosion energy of the supernovae based 
upon nucleosynthesis, measured neutron star masses, and light curves are discussed 
and applied. The nucleosynthetic results, when integrated over a Salpeter initial 
mass function (IMF), agree quite well with what is seen in the sun. We discuss in 
some detail the production of the long lived radioactivities, ^^Al and ^'^Fe, and why 
recent model-based estimates of the ratio ^^Fe/'^^Al are overly large compared with 
what satellites have observed. A major source of the discrepancy is the uncertain 
nuclear cross sections for the creation and destruction of these unstable isotopes. 



1 Introduction 



Starting in the late 1970's, with the encouragement of his good friend Gerry 
Brown, Hans Bethe became interested in applying his expertise in nuclear 
physics to one of the more vexing problems in modern astrophysics - how 
massive stars die as supernovae. The problem is difficult for a variety of rea- 
sons. The iron core of a massive star collapses to a neutron star (or sometimes 
a black hole) and, somehow, some fraction of that remnant's binding energy 
is converted into outwards kinetic energy in the overlying star. The favored 
model, now as then, says that the binding energy of the neutron star is radi- 
ated as neutrinos, a fraction of which deposit their energy in the matter above 
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the neutron star causing it to expand and explode (jColgate fc White 



19661). 



When Hans began to work on the problem, supernova models were not giving 
explosions. Moreover, the physics was very unc ertain w i th bo unce densities 



ranging from 10^^ g cm""^ to 10^^ g cm~^ (e.g., IWilsoru . Il978l ). The nuclear 
equation of state was partic ularly uncertain. A major breakthrough was the 
work by iBethe et al.l (119791 ) who showed that the heat capaci ty of the nu- 



clear bound states was much larger than previously believed (IFowler et al. 
19781 ). In fact, the mean excitation energy was ~ {A/8){kT)'^ and the 



partition function associated with all these states was exponentially huge, 
G ~ exp((A/8) (kT)). Consequently, nuclear equilibrium favored bound nuclei 
which remained abundant, increasing their average mass, until they touched 
and merged near nuclear density. Bounce occurred at super-nuclear density on 
the hard core, repulsive component of the strong force (not thermal pressure 
as some calculations claimed) and was at low entropy. The general idea of 
entropy as an important variable in core collapse came from Hans, who liked 
to remark that though the bounce was thermally very hot, in terms of entropy 
it was as ordered as ice. 



During the next 20 years, Hans made many other lasting contributions to 
supernova theory, including the currently favored "de layed" neutrino t r anspo rt 
paradigm in which convection plays a major role teethe fc Wilson! . ^dMT 
Hans also introduced the ideas of a "gain radius" , where neutrino heating first 
exceeds neutrino losses, and of "net ram" , the momentum of the accretion flux 
that must be overcome to get the shock to move out. He excelled in simple 
analytic models for the physics of core collapse, and brought a much needed 
physically intui tive understa nding of a subject that had hitherto been largely 



numerical (e.g., Bethe, 1990l ). 



Because of his historical interests in stellar structure and nuclear physics, 
Hans was also interested in the presupernova evolution of massive stars and 
in nucleosynthesis. During his visits to Santa Cruz and by mail, we had many 
discussions on the progenitor of SN 1987A, the physics of supernova light 
curves, the nature of the "reverse shock", explosive nucleosynthesis, and on 
the r-process. Thus it is to his memory that this paper is dedicated. 



To th is day, we still don't know exactly how massive stars explode flWoosley fc Janka 
20051 ). so the parameterization of the explosion is discussed in §3. The key nu- 



^ The first calculation to show the revival of the shock by neutrino heating was 
carried out by Wilson alone in 1982 (jWilsonl . Il985l ). but analysis of the calculation 



and the first refereed publication was by B ethe and Wi l son. F or a time, Hans also 
embraced the idea of "prompt explosions" (jBaron et al.l . 119871 ) , explosions in which 
neutrino transport played no constructive role and the explosion was due to a hy- 
drodynamical "bounce" . He gave up the idea after detailed calculations showed that 
neutrino losses and photodisintegration killed the prompt shock. 
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clear reaction rates and other uncertain aspects of the presupernova evolution 
are described in §2. In §4, our principal nucleosynthetic results are presented, 
and in §5, we conclude with a discussion of two key species of interest in 7-line 
astronomy, ^^Al and ^°Fe. 

Throughout this paper and in the future, we employ a unit of energy, the 
"Bethe", abbreviated "B", equal to 1.0 x 10^^ erg. Gerry Brown introduced, 
and Hans and Gerry both promoted the use of an alternate term, "foe", fre- 
quently found in the supernova literature to stand for 10^^ erg, but in defer- 
en ce to Hans ' cont ributions to the field, we follow the convention suggested 



bylWeinbergl (120061 ) 



2 Uncertainties in the Presupernova Evolution 



2. 1 Critical Reaction Rates 



The key uncertain reaction rate affecting both the structure of and nucle- 
osynthesis in massive stars remains ^^Cf a, 7)^^0, despite over 30 years of 



painstaking laboratory investigation fe.g.. iDyer &: Barnes. Il974f) . The exper- 
imental situation was recently reviewed by iBuchmanru (120051 ). who recom- 
mends S(300 keV) = 102 - 198 keV b with a best value of 145 keVb. Based 
upon nucleosynthesis considerations, IWeaver fc Woosleyl (119931 ) estimated an 
S-factor of ~170 keVb, which remains within the experimental range today. 
More precisely, Weayer an d Woosley suggested a rate 1.7 ± 0.5 times that of 



Caughlan fc Fowled (119881 ). which would be 120 - 220 k eVb, but even at the 



time, the error bar was regarded as liberal. More recently. iBoyse. Woosley. Sz Heger 
( 120021 ) has revised the nucleosynthesis constraints using more stellar models, 
a finer grid of ^^C{a, 7)^^0 rates, finer stellar zoning, and other improvements 
to the stellar model. Their results, shown in Fig. 1, are in good agreement 
with the earlier calculations of Weaver and Woosley, but give a narrower error 
bar and also make the sensitivity of the results to this rate (variations of only 
10% matter) more apparent. Because of the need to include a rate that is 
accurate across a wide range of temperature, not just during helium burning, 
the preferred rate is again ex pressed as a multiple of a published rate fit, this 



time 



BuchmannI (119961 . 119971 ). which has S(300 keV) = 146 keVb. Boyes' best 
fit is about 1.2 times this, or 175 keVb and a value of 1.2 times Buchmann 
(1996) was used in the present study. This is also consistent with recent mea- 
surements reported by iHammer et al. (l2005h that give a best value of 1.08 
times Buchmann-1996 (i.e., 162 ±39 keVb). Of course, one could argue that 
the nucleos ynthesis limit is also influe nced by our uncertain model of stellar 



convection (jWeaver &: Woosleyl . Il993l ). in which case an experimental value 



ultimately near 170 keVb would serve to validate the treatment of convection 
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Fig. 1. Production factor factor of key elements for a set of solar met allicity stars 
folded with a Saltpeter birth function (jBovse. Wooslev. &: Hegerl . |2002| ) . 



in the code. 



During the end of helium burning the ^^C(a,7)^^0 rate competes with the 
triple-alpha reaction rate, and hence the uncertainty in that rate can have 
similar effects. In test calculations at 3 x 10^ K and 1000 and 2000 g/cm^ 
we found that a 10 % increase in the triple-alpha rate would have the same 
consequence as an 8% decrease in ^^C(a;, 7)^^0. In a star, convection may 
change these results, though probably not much. The ^^C(a,7)^^0 rate would 
need to be known better than abou t 10 % before the uncertainty in the triple 
alpha rate, ~12 % flTur et al.l . 120071 ). becomes a limiting factor. 



The other uncertain reaction rate that affects the abundances of hosts of nuclei, 
not just a few, is ^^Ne(a,n)^^Mg, which, along with ^^Mg(n,7)^^Mg, regulates 
the strength of the s-proce ss in massive stars. The rate used here is the rec- 
ommended value from from I Jaeger et al.l (120011 ). The reaction ^^Ne(a, 7)^^Mg 
competes with ^^Ne(Q!,n)^^Mg for the destruction of ^^Ne and is thus also of 
some importance for determining the strength of the s-process. Here we use 
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the lower bound for ^^Ne(a, 7)^^Mg estimated by lKappeler et al.l (Il994r). Other 



choice s of strong and weak reaction rates have been discussed by lWoosley et al. 



(120021 ). In particular, except where o t herwi se noted, we use the Hauser-Feshbach 
rates from Rauscher fc Thielemann fl2000h for reactions involving n, p, and a 
on heavy nuclei that lack experimental determination. This is of some rele- 
vance to the issue of ^^Al and ^°Fe production discussed later in the paper 
(^5.1). 



2.2 Mass Loss 



Mass loss is known to be a powerful determinant in the evolution of stars of 
nearly solar metallicity, and its om ission was one of the maj or shortcomings 
of the Woosley- Weaver 1995 survey IWoosley fc Weaver! fjl995l ). For stars more 
massive than about 35 Mq, mass loss is particularly important since it not 
only removes the hydrogen envelope, but shrinks the helium core appreciably. 
With current estimates of mass loss, a 100 Mq Population I star ends its life as 
a star of only about 6 Mq, composed of helium and heavier elements only and 
no hydrogen left. This is similar to the mass and composition of the core of a 
20 Mq star, and the explosion properties, remnant mass, and nucleosynthesis 
are radically different from a 100 Mq star that had no mass loss (see §3.1). 



The m ass loss prescription us ed here has also been discussed by Woosley et al. 



(I2OO2I ). In particular, we use iNieuwenhuijzen fc DeJagerl (119901) fo r mass loss 
on the main sequence and for red giants and Wellsein fc Langerl (19991) for 
Wolf-Rayet stars The latter is based on the mass loss rate by iBrauru (119971 ) fit 
to observational data a nd div ided by a factor of three to account for clumping 
Hamann fc Koesterkel (Il998l ). The nucleosynthesis products carried away by 
stellar wind are included in all yields reported in this paper. 



2.3 Convection and Rotation 



The treatment of convective phy sics, including overshoot m ixin g and semicon 



vection, follows the discussion in I Woosley Sz Weaver! (119881 ) and lWoosley et al 



(120021 ). In particular, we use a semi-convective mixing parameter, a = 0.1, 
which results in relatively fast mixing in semiconvective regions. Mixing was 
treated in a time-dependent, mixing length formalism using the Ledoux crite- 
rion for instability. The fast semiconvection contributes significantly to mix- 
ing in regions that are stable to the Ledoux criterion but unstable to the 
Schwarzschild criterion, however, the mixing is less than that of a mere Schwarzschild 
mixing, taking into account the stabilizing effects of composition gradients. 



Rotation can have a significant effect on both the presupernova evolution and 
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the explosion mechanism. Here rotation was neglected, which is to say it is 
assumed that the change in helium core mass and dredge up of light isotopes 
due to rotationally induced mixing are small, for moderately fast rotating 
stars, and that the ratio of centrifugal force to gravity during the explosion is 
negligible. All these assumptions are questionable for rapidly rotating stars, 
espec ially so in small fraction of massive stars that become gamma-ray bursts 
(e.g., IWoosley fc Hegerl . I2OO6I ). 



2.4 Initial Abundances 



The evolution and nucleosynthesis of a massive star both are sensitive to its 
initial composition. The total abundance of CNO affects the efficiency of hy- 
drogen burning and the opacity. The conversion of CNO into ^^Ne during 
helium burning determines the "neutron excess", which affects the production 
of all nuclei with unequal neutron and proton numbers. ^^Ne also provides the 
free neutrons necessary for the s-process during helium burning. Finally, be- 
cause the yields of supernovae are traditionally normalized to ^^O, any change 
in the solar oxygen abundance affects the comparative ease with which all 
other heavy elements are produced. 



It is thus a major occurrence in nucleosynthesis theory when the solar abun- 
dances, traditionally taken as representative of Population I stars in our Galaxy, 
ar e modified. Rece nt re visions to the solar abundance set h ave been discussed 
by iLoddersI (120031 ) and lAsplund. Grevesse. fc Sauvall (120041 ) . The abundances 
of all isotopes of CNO have been reduced by am ounts of order 30 % compa red 
with the standard Anders and G revesse v a lues (Anders fc Grevessd . Il989l ) of 
a few years ago. Here we use the iLoddersI (120031 ) set both as a starting com- 
position, and also to normalize all computed yields. 



2.5 Presupernova Models 



Using the Kepler implicit h ydrodynamics code ([Weaver et al.l . Il978l ) and the 
physics specified above and in lWoosley et al.l (120021 ). stars of solar composition 
and various masses were evolved to the presupernova stage - defined by a 
collapse velocity in the core of 1000 kms~^. Masses included in the study were 
12 through 33 solar masses in steps of 1 Mq, plus stars of 35, 40, 45, 50, 55, 60, 
70, 80, 100, and 120 Mq - 32 stars altogether. A future survey will use a much 
finer grid of masses, and the present work may be regarded as a preliminary 
survey. 
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3 Simulating the Explosion 



As alluded to in the introduction, a robust description for how massive stars 
explode as supernovae remains elusive and this must surely affect our under- 
standing of the origin of the elements. It is worth separating out that part of 
the nucleosynthesis that depends on the explosion mechanism from that which 
does not, however. 

Certainly isotopes produced in the vicinity of what is commonly known as "the 
mass cut" are sensitive to conditions set up by the passage of the shock. This 
includes the yields of species made by explosive oxygen and silicon burning and 
in nuclear statistical equilibrium. More quantitatively, these are the species 
made at temperatures above 3 x 10^ K and at radii less than about 7000 km, 
i.e., roughly the inner 1 to 2 solar masses of ejecta. Other species made by 
hydrogen, helium, carbon, neon and oxygen burning in hydrostatic equilibrium 
are not greatly affected (provided such material escapes the star and does not 
fall into a black hole), nor is the nucleosynthesis in the pre-explosive wind. On 
the other hand, the r-process and other species made in the neutrino-powered 
wind are quite sensitive to the explosion mechanism, and it is this sensitivity 
that makes them excellent diagnostics of the event. 

As we shall see though, even the "explosive nucleosynthesis" below atomic 
mass 100 is not particularly sensitive to details of the explosion, provided that 
the star blows up with a "reasonable" kinetic energy and the explosion is not 
grossly asymmetric. This is basically because the shock conditions are deter- 
mined by the pre-explosive structure and some simple physics, AirR^aT^/S = 
explosion energy ^ 1 B. 



Here, as elsewhere, the explosion is parameterized by a piston at constant 
Lagrangian m ass coordinate that moves through the star with some specified 
radial history (IWoosley fc Weaverl . Il995l : IWoosley et al.l . |2002| ) . The essential 
parameters of the piston are its location in mass and the final kinetic energy it 
imparts to the ejecta at infinity. Two different choices of each are explored: a) 
piston mass at the edge of the iron core or at the point where the dimensionless 
entropy S/N^k = 4.0; and b) kinetic energies of 1.2 and 2.4 B. Thus for 
each mass, 4 explosion models were calculated for a total of 128 supernovae 
simulated. 



The choices of piston mass and explosion energy are not free parameters, but 
are highly constrained by observations. The piston mass cannot be smaller 
than the iron core mass or unacceptable overproductions of ^^'^^Fe and other 
neutron-rich species in the iron group will occur. On the other hand it cannot 
be much larger than the base of the oxygen shell {S/NAk = 4) or, as we shall 
see, typical neutron star masses will be too large. The large density decrease 
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Fig. 2. Entropy and density distributions inside a 20 solar mass presupernova star. 
The iron core mass here is 1.54 Mq; the base of the oxygen shell is at 1.82 Mq. The 
sudden decrease in density at the base of the oxygen shell causes an abrupt decline 
in ram pressure which often results in explosions happening with this mass cut. 

associated with the base of the oxygen shell is also dynamically important and 
successful explosion calculations, when they occur, frequently find the mass 
cut there. The explosion energy is constrain ed to be 1 - 2 B by observations of 
SN 1987A (lBethel . ll99nl : lArnett et all . ll989l ) which was a Type II supernova of 
typical mass (about 18 - 20 Mq). It is also constrained by the observed light 
curves of Type II supernovae. 



3. 1 Remnant Masses 



Observations by iThorsett fc Chakrabartyl (119991 ) of a large number of pulsars 
in binary systems give a narrow spread in masses, 1.35 =fc 0.0 4 Mq. There 
must be room for some diversity, however. iRansom et al.l (120051 ) present com- 
pelling evidence for a pulsar in the Terzian 5 globular cluster with a gravi- 
tational mass of 1.68 Mq. The remnant gravitational masses for our survey 
using the Kepler stellar evolution code, with KE = 1.2 B and pistons lo- 
cated the edge of the iron core, are plotted in Fig. 3. A more careful analysis 
of fall back in these models using an Eulerian hydrodynamics code and a 
better treatment of the inner b oundry conditions has been carried out by 



Zhang. Woosley. fc Hegerl (120071), but gives simila r numbers for solar metal- 



licity stars. Using the IZhang. Woosley. &: Hegerl (120071 ) values, adopting a 
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Fig. 3. Fe core masses for a grid of stellar masses. Solar metallicity stars. See text 
for explanation. 

Salpeter initial mass function with F = 1.35 to describe the birth frequency 
of these stars, and assuming a maximum neutron star mass of 2.0 Mq, one 
obtains an average gravitational mass for the neutron star of 1.47 ± 0.21 if 
the piston is at the S/N^k = 4.0 point and 1.40 ± 0.22 if it is at the edge 
of the iron core. If the maximum neutron star mass is 1.7 Mq, the numbers 
are changed to 1.41 ± 0.15 Mq and 1.34 ± 0.14 respectively. In this paper, we 
carried out simulations with the piston at both points - the iron core edge, and 
the base of the oxygen shell. Larger masses than typical are also possible for 
the rare exceptionally massive star, usually those over 25 Mq. For those cases 
where a black hole was made, its average mass was around 3 Mq. We note 
that these numbers are for single stars and they could be altered significantly 
in mass exchanging binaries. 



The figure also shows that neutron stars are made by both the lightest main 
sequence stars and the heaviest. This is a consequence of mass loss. The helium 
core mass of the presupernova star increases monotonically with main sequence 
mass up to about 45 Mq, where it reaches a maximum of 13 Mq. Beyond that 
the helium core shrinks due to efficient Wolf-Rayet mass loss and the iron core 
mass shrinks with it. A 100 Mq model had a total mass of only 6.04 Mq when 
it died - all helium and heavy elements - and an iron core mass of 1.54 Mq. 

The results are quite different for stars with low metallicity and, hence, reduced 
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Fig. 4. Fe core masses for a grid of stellar masses. Zero metallicity stars. Many more 
black holes are made because the star loses little mass during its evolution to pre- 
supernova. The two branches of black hole masses at high main sequence mass cor- 
respond to red (lower branch) and blue (upper branch) supergiant progenitors. The 
lower carbon-oxygen core mass es for the red supergiant c a ses re f lect dredge up and 
prini ary nitrogen production (IZhang. Woosley. h Hegerl . 120071 ; iHeger &: Woosleyl . 
200?! ^. 



mass loss (IHeger fc Woosleyl . 120071 : IZhang. Woosley. fc Hegerl . l2007l ). Fig. 4 
shows that the remnant mass increases rapidly for main sequence masses above 
about 35 Mq and continues to increase at higher masses. These large masses 
are due to fall back. A 1.2 B explosion is inadequate to unbi nd the entire star, 
especially given the large helium core (IWoosley et al.l . |2002|) and effect of the 
reverse shock. A 100 Mq main sequence star now dies wit h a helium core of 
42 M0, well into the pulsational pair instability domain (IHeger &: Woosleyl . 



2OO2I ). Unless supernova engines of much greater power than 1.2 B become 
available at low metallicity, these stars will make black holes, not neutron 
stars, and if the rotation rate is sufficient, gamma-ray bursts. 

One may also not e the existence of two branches of black hole remnants above 
35 Mq in Fig. 4. IZhang. Woosley. fc Hegerl (120071 ) find that these branches 
correspond to two different classes of progenitors - red supergiari t s, wh ich ex- 
perience a lot less fall back during the reverse shock (IChevalierl . Il989l ) - and 
more compact, extremely blue supergiants. If the star produces primary nitro- 
gen due to the interpenetration of the helium convective core and hydrogen 
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envelope, it swells to red giant proportions, has a weaker reverse shock, and 
leaves a smaller remnant mass. 



3.2 Light Curves 



The KEPLER code includes radiative diffusion and can thus be used to cal- 
culate approximate light curves for the supernovae it produces. The code is 
limited by using a single temperature for the radiation and the matter, and 
assumes blackbody radiation, but th ese are not bad approximations during the 



plate au stage of Type 11 supernovae (jWeaver fc Woosleyl . Il980l : Eastman et al. 



19941 ). The principal opacity source is electron scattering with the free electron 



abundance determined by solving the Saha abundanc e s of a ll ions for the 19 
isotopes in the reaction network (lEnsman fc Woosleyl . Il988l ). A floor opacity 
of 10~^ cm^ g~^ is used in regions that have recombined. The abundance of 
^^Ni is taken from the nucleosynthesis calculation and moderate mixing of the 
helium core is assumed. 



The resulting light curves for four explosions of a 15 M© supernova are given in 
Fig. 5 for cases where the mass cut was taken at the edge of the iron core and 
at the location where the entropy equals 4/cB/baryon. Two explosion energies, 
1.2 B and 2.4 B were employed. The explosions that had the higher kinetic 
energy were brighter on the plateau and the ones with the deeper mass cut 
(and hence more ^^Ni ejected) had the brighter tails. The mass of ^^Ni ejected 
was 0.086 Mq for the 1.2 B explosion with the mass cut at S = 4/cB/baryon; 
0.096 M0for the 2.4 B explosion with mass cut at S = 4A;B/baryon; 0.27 Mq 
for the 1.2 B explosion with the mass cut at the edge of the iron core; and 0.31 
Mq for the 2.4 B explosion with the mass cut at the edge of the iron core. 



Cle arly, the mode ls with higher kinetic energy are brighter on the plateau (see 
also iPopovl . Il993l ) . In fact, if the kinetic energy were any larger than 2.4 B, the 
supernova would be far brighter than average Type lip supernovae. On the 
other hand if the explosion energy were much less than 1 B, large amounts of 
material would fall back, increasing the masses of the neutron star remnants 
beyond acceptable values and robbing the nucleosynthesis of its most prolific 
sources. We conclude that the range 1.2 - 2.4 B is the relevant one for modern 
supernovae in solar metallicity stars and these are the values employed in the 
nucleosynthesis survey. 
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Fig. 5. Light curve of a 15 solar mass supernova. The four curves represent two 
choices of explosion energy and combined with two choices of mass cut. The presu- 
pernova star was a red supergiant. 

4 Nucleosynthetic Yields 



The integrated yields of the elements are given for four different choices of 
mass cut and explosion energy in Fig. 6. Whether one places the piston at the 
edge of the iron core or the base of the oxygen shell and whether the explosion 
energy is 1.2 B or 2.4 B makes little difference except to the iron group. There 
the difference is of order a factor of two, with lower iron yields obviously 
resulting from lower explosion energies and shallower pistons. In all cases the 
iron group synthesis is low compared both with C, O, Ne, and Na and with 
s-process production above Ni. One expects fro m one-half to t wo-th irds of 



the iron group to come from Type la supernovae (ITimmes et al.l . Il995l ) which 
are not included here. The s-process, which is secondary in nature, will be 
underproduced in stars of less than solar metallicity, so a factor of two extra 
here relative to oxygen is not undesirable. 

Fig. 7 shows the integrated nucleosynthesis (the yields folded with a Salpeter 
initial mass function with F = -1.35) for all elements up to lead compared with 
the isotopic composition of the sun. Fig. 8 shows the corresponding comparison 
of isotopes. Overall, the agreement is quite good, especially considering that 
several known sites of important nucleosynthesis have been omitted. Classical 
novae will need to produce ^^N and ^^O, though some ^^N is made here by 
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Fig. 6. Elemental yields integrated over a Salpeter initial mass function for solar 
metallicity stars with masses from 12 Mq to 120 Mq. Calculations were carried out 
for for two choices of explosion energy (1.2, 2.4 B) and two piston location (at the 
place where the entropy jumps to S/NAk = 4 and at the edge of the iron core). Only 
minor differences are discernible. Similarly small changes occur when the slope of 
the IMF is changed from -1.35 (Salpeter) to a gentler -0.9. 

the neutrino spallation of ^^O. The isotopes ^^Ca, '^^Ti, a nd ^^Ca are und er- 
produced and may come from some rare form of SN la (jWoosIeyl . 119971 ) or 
asymmetric supernova. The overabundance of is not a concern since some 
will decay before the sun is born. Also missing are ordinary Type la super- 
novae, which contribute half or more of the iron group, and asymptotic giant 
branch stars which make ^^N and the s-process. In fact, the full production 
of carbon here is a novel and surprising result, since it is usually attributed 
to low mass stars. It is made here chiefly in the winds of very massive Wolf- 
Rayet stars and its production is facilitated by the new lower solar abundance 
(jLoddersl . [2nn3l ). 



Finally conspicuously absent is the r-process and other products of the neutrino- 
powered wind. The wind of a young neutron star is a prolific source of heavy 
elements, accounting fo r about half of the is otopes in nature. These include 



not onl y the r-process (IWooslev et al. 



nuclei (iPruet et al.l. 12006 : iFrohlich et al 



like Zn ( iHoffman et al 



19961 ) and Sc ( 



■994r). but some important p-process 



20061). and even abundant elements 



Pruet et al.l . |2005| ). Hans Bethe was 



very interested in the neutrino-powered wind and the r-process, and he was 
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20 40 60 

element charge number 



Fig. 7. Elemental yields when the grid of supernova masses is integrated over a 
Salpeter initial mass function. The explosion energy was 1.2 B and the piston was 
located at the mass coordinate where S/iV^/c = 4.0. A strong s-process operates up 
to Z = 40. 

working on it when SEW last saw him in Winter 2003. This was probably 
his last supernova-related project. He said that he had an abiding interest in 
uranium. 



5 The Special Cases of ^^Al and ^^Fe 



Having computed the isotopic nucleosynthesis in a grid of stars up to 120 
Mq, including the contribution of the winds the more massive stars make 
as Wolf-Rayet stars, we turn to the examination of two isotopes of special 
interest to gamma-ray line astronomy. The long-lived isotopes ^^Al and ^°Fe 
accumulate in the interstellar medium from thousands of supernovae and thus 
serve as cali bration s on t he integrated yields of massive s tars. Observations 
by RHESSI fISmithl . \2004 \ and INTEGRAL flHarris et all . l2005h give a ratio 
of fluxes from the decays of ^^Al and ^°Fe of 0.16 and 0.11 ± 0.03 respectively. 
Both measurements are quite consistent with th e predicted value, . 16, by 



Timmes et al.l (119951 ) based upon yields from the IWooslev fc Weaverl (119951) 



survey. Later calculations (iRauscher et al.l . 120021 : iLimongi fc Chieffil . l2003h 
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Fig. 8. Isotopic nucleosynthesis integrated over a Salpeter initial mass function. 
The results are in good agreement with solar abundances below A = 85. The excess 
of s-process nuclei above the iron group is needed in part, to compensate for the 
smaller production of these secondary nuclei in stars of slightly lower metallicity. 
The 7-process is mostly successful in making the neutron-deficient isotopes (the 
"p-nuclei") above A = 130, but there is an annoying deficiency of p-nuclei production 
from A = 90 to 130. The large productions of ^^B, ^^F, ^^^La, and ^^^Ta, are due to 
the neutrino process. Products of the neutrino wind, e.g., ^^Zn and the r-process, 
are not included in this plot. 
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however, using stellar and nuclear physics that was nominally "improved" 
gave a much lar ger synthesis of ^''Fe that was not in hne with observations 
flPrantzosl . [2003 ). 



The ratio of gamma-hne fluxes implies, in steady state, a synthesis ratio by 
mass of ^°Fe/^^Al of 60/26 times the flux ratio, or about 0.3. (The steady 
state abundance is inversely proportional to the decay rate and the flux is the 
abundance times the decay rate so the decay rate itself cancels.) iTimmes et al. 
( I1995I ) gave a theoretical ratio of 0.38 with an expected uncertainty of a factor 
of 1.7. Using the larger grid of models here, however, and including mass loss 
as discussed in §2.2, we calculate a ratio of 1.8, i.e., six times t o o larg e. This 
large excess of ®°Fe/^^Al is consistent with what iRauscher et al.l (120021 ) found, 
even though their calculations did not include the quite massive stars studied 
here (up to 120 Mq), nor their mass loss. What is wrong? 



5.1 Nuclear Physics Uncertainties 



One problem is certainly the use of uncertain nuclear reaction rates in all 
studies to date. In mak i ng the transition to the reaction rate data base of 
Rauscher fc Thielemann! (120001 ). we erroneously used the new Hauser Fes- 
hbach rates especially for ^^Al(n,p)^^Mg and ^^Al(n,«)^^Na. These are the 
principal means of ^^Al destruction in the carbon and neon layers where ^^Al 
is explosively synthesized. The Rauscher- Thielemann rate for ^^Al(n,p)^^Al 
at 2 X 10^ K, for example, is 1.4 x 10® cm^ Mole^^ s"^. The Rauscher- 
Thielemann rate for '^^Al{n,ay^Na at 2 x 10^ K is 2.6 x 10^ cm^ Mole"^ 
s~^. These a re a both a factor o f 3 to 5 higher than the rates used for these 
re actions bv IWooslev fc Weayerl (I19951) and the experim ental determinations 
by iKoehler et al.l (Il997h and (ICaughlan fc Fowled Il988h . 



A second effect, less important than the cross sections, is a super-hot helium 
shell (4 x 10® K) in several of the pre-supernova star. This shells existence 
was traced to the use of OPAL opacities in a region where they may not be 
appropriate, a region where electron scattering dominates. Using the electron 
scattering opacity of I Weaver et al.l (119781 ) just in high temperature regions 
where electron scattering dominates decreased the ^°Fe yield significantly, but 
this was only in a few stars. 

Using what we believe to be more nearly correct cross sections for ^^Al de- 
struction (though still uncertain) and adjusting the opacity as described, the 
integrated yield of ^°Fe to ^^Al is reduced to 0.95. This is for a standard 
Salpeter IMF with F = -1.35. If we instead change the slope to -0.90, i.e., 
enhance the production of very massive stars, the ratio is reduced slightly to 
0.81. Even then one-half the yield of ^^Al comes from stars under 35 Mq, not 
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the more massive ones and their Wolf-Rayet winds. 



There are further uncertainties to explore, however. Several of the cross sec- 
tions governing the production of ^°Fe are also highly uncertain. The reaction 
^^Fe(n,7)^°Fe affects its synthesis and ^''Fe(n,7)^^Fe controls its destruction. 
Neither is measured, thoug h both could be, admitted l y with difficulty. Interest- 
ingly, both changed in the Rauscher fc Thielemann ( 2000l ) tab ulation in such 



a dire ction as to increase Fe production. The tabulation by IWoosley et al. 



(119781 ) had, for helium burning temperatures, a rate for ^^Fe(n,7)^°Fe half as 
large and a rate for ®°Fe(n,7)^^Fe twice as large. When the older rates were 
used for a select set of models, the ®°Fe production was reduced by about a 
factor of two. 

The final nuclear uncertainty is the rate governing the production of neutrons 
where ^°Fe is made , i.e., ^^Ne(a,n)^^Mg. The rate included in our network 



(jjaeger et al.l . 1200 ll ) is increased from what was used in 1995. If we reduce its 
value in a few select models by a factor of two (within the error bar), ^°Fe 
production is again decreased by up to a factor of two, though usually the 
effect was smaller. 

All things considered, variation of only the nuclear physics, bringing uncertain 
cross sections back to the values they had in the Timmes et al survey, could 
account for most of the difference in the present calculations and the obser- 
vations. Hence further progress in this important field of astronomy depends 
upon more accurate measurements and estimates of critical nuclear physics. 



5.2 Uncertainties in the Stellar Models 



This is not to say that non-nuclear effects are unimportant. Met allicity, mass 



l oss, r otation, and an uncertain IMF certainly all play major roles. iPalacios et al. 
(120051 ) have explored ^^Al production in models of massive stars that include 
rotationally induced mixing, as well as mass loss and different choices of metal- 
licity. An explicit comparison with a couple of our models is educational. For 
a 60 Mq main sequence star with Z = 0.02 and no rotation, they find an ^^Al 
production in the pre-explosive wind of the star of 1.30 x 10~^ and a final star 
mass of 12.4 M© (iMevnet fc Maederl . [200o[ ) . For our 60 Mq model with metal- 
licity Z = 0.016 and using the smaller experimental rates for ^^Al(n,p)^^Mg 
and ^^Al(n,a)^^Na, we find a production in the wind of 1.1 x 10~^ M, 
final mass of 8.0 Mq. But we also find an additional 9.9 x 10"^ Mq 
is produced in the explosion of the remaining star, chiefly by explosive neon 
burning. This is good agreement, and shows that the explosion and wind may 
contribute cornparab le amounts to ^^Al synthesis even for a 60 Mq progenitor. 



Q and a 

of 26 Al 



Palacios et al.l (120051 1 further explore the dependence of metallicity and rota- 
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tion though, and find ^^Al production in the wind of this same star is increased 
to 2.2 X 10~^ Mq if the rotation rate is 300 km s~^ on the main sequence, or 
3.0 X 10~^ Mq with no rotation but Z = 0.04. Combining both effects, Z = 
0.04 and v^ot = 300 km s~^, the ^^Al production in the wind becomes even 
larger 7.2 x 10"^ Mq. While one must be concerned that increasing the metal- 
l icity may also in crease the ^°Fe yield and thus not increase the ^°Fe/^^Al ratio 
(jPrantzosl . |2004| ) . this does show the sensitivity of ^^Al to reasonable variations 
in rotation rate. For a 120 Mf^ r nodel , the effect is even greater is similar. For 
Z = 0, v,ot = 0, IPalacios et all fl2005h obtain an ^^Al mass of 5.7 x 10"^ Mq 



in the wind while we have 4.9 x 10 ^ Mq plus 2 .9 x 10 ^ M^d made in the 
explosion. With Z = 0.04 and V^ot = 300 km s"\ IPalacios et aU (j2005l ) get a 
whopping 2.2 X 10^3 Mq. 



Limongi Sz Chieffil (120061 ) have also recently (after our present study was com- 
pleted) examined the sensitivity of ®°Fe and ^^Al production to the prescrip- 
tion for mass loss and the slope of the IMF. They find that both can make an 
appreciable difference. 



6 Conclusions 



We still don't understand exactly how massive stars explode, far less the vari- 
ation of explosion properties - especially mass cut and explosive kinetic energy 
- with main sequence mass. This remains a forefront problem in nuclear as- 
trophysics research to which Hans Bethe contributed greatly. It is likely, in 
the final analysis, that the physical intuition, terminology, and convective, 
neutrino-powered paradigm that he and his colleagues brought to the field 
will form the basis of a c omplete understanding, though we aren't there yet 
(IWoosley fc Jankal . l2005l ). Certainly, the low entropy, super-nuclear density 
bounce following the initial collapse will be phase one of any massive star 
explosion. 



This lack of a first principles understanding of the explosion mechanism, how- 
ever, is not a fundamental roadblock on our path to understanding the origin 
of (almost all of) the elements. Arguments have been presented here to show 
that the mass cut is highly constrained by nucleosynthesis and observed neu- 
tron star masses. The explosion energy in common Type 11 supernovae is also 
mostly in the range 1.2 B plus or minus a factor of two. Exploding a large 
range of stellar masses with pistons located either at the edge of the iron core 
or the base of the oxygen burning shell - the maximum range allowed - and 
with explosion energies of either 1.2 B or 2.4 B gives very similar nucleosyn- 
thesis. The iron group is most affected and the magnitude of the uncertainty 
is about a factor of two. 



18 



The nucleosynthesis that results (Fig. 7 and Fig. 8) agrees reasonably well with 
solar abundances. There are some changes caused by the recent downward 
revisions of the solar CNO abundances, and at first glance the agreement is 
worsened by these changes. Since ^^O is our standard normalization point in 
nucleosynthesis studies, since we now need to make less of it, the production 
of all other heavier elements is decreased. Yields that previously would have 
coproduced Si and O say, in solar proportions, now overproduce O (Fig. 8). The 
production of odd-Z elements and odd-A isotopes is also decreased because 
the initial CNO in the star later becomes the ^^Ne that sets the neutron excess 
for carbon, neon, and oxygen burning (§2.4). Still the agreement is not too 
bad, and most of the missing species ^^C, ^^'^^N, ^^Ca, etc. can be attributed 
to other sites than massive stars. 



The outstanding problem in nucleosynthesis theory presently is a full under- 
standing of the r- and p— processes. The latter has an appreciable contribution 
from explosive neon and oxygen burning (shown in Fig. 8) for A greater than 
130, but is underproduced for lighter masses. The solution for both the r- 
process and the light p-process probably lies in the neutrino-powered wind. 
Current models give inadequate entropy in the wind and this may be where 
nucleosynthesi s can be an important diagnostic of the explosion model (e.g.. 



Burrows et al.l . 120061 ). 



The nucleosynthesis of the long-lived radioactivities ^^Al and ^°Fe is an im- 
portant constraint on the stellar models, and one that is largely independent 
of the explosion mechanism. The abundances inferred from gamma-ray line 
astronomy may have important implications for rotationally induced mixing, 
convection theory, mass loss theory, the initial mass function for massive stars, 
and the distribution of metals in the galaxy. The synthesis is also quite sen- 
sitive to nuclear reaction rates whose uncertain values could be better de- 
termined in the laboratory, however. In particular, the discrepancy between 
o bservations of the ^° Fe/^^Al ratio and recent calculations - this work and 



Rauscher et al.l (120021 ) - may involve a "perfect nuclear storm" of erroneous 
choices. The rates affecting ^^Al destruction were almost certainly too high; 
the rates affecting ^°Fe production, namely ^^Fe(n,7)®°Fe and ^^Ne(a,n)^^Mg, 
may have been too high; and the rate for its destruction, ^°Fe(n,7)^^Fe may 
have been to o low. Given the choi ces made by Woosley fc Weaver! ( 19951 ). the 
prediction of iTimmes et al.l (119951 ). which agrees with observations, is still de- 
fensible. In any case, important inferences about the stellar models will only 
be credible (and necessary), when these uncertain rates have been better de- 
termined. 
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